Electron-phonon coupling in the Cqq fullerene within the many-body GW approach 
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We study the electron-phonon coupling in the 6*60 fullerene within the first-principles GW ap- 
proach, focusing on the lowest unoccupied ti u three-fold electronic state which is relevant for the 
superconducting transition in electron doped fullerides. It is shown that the strength of the cou- 
pling is significantly enhanced as compared to standard density functional theory calculations with 
(semi)local functionals, with a 48% increase of the electron-phonon potential V ep with respect to 
the LDA value. The calculated GW value for the contribution from the H g modes of 93 meV comes 
within 4% of the most recent experimental values. The present results call for a reinvestigation of 
previous density functional based calculations of electron-phonon coupling in covalent systems in 
general. 

PACS numbers: 71.15.Qc, 71.20.Tx, 74.20.Fg, 74.20.Pq 



I. INTRODUCTION 

Electron-phonon coupling in molecular systems is 
at the heart of several important physical phenom- 
ena, including the mobility of carriers in organic elec- 
tronic devices? 1 - - — the dissociation of excitons at the 
donor/acceptor interface in organic photovoltaic cells,— 
or the superconducting transition in molecular solids, 
from the most famous fulleride case£~— to the recent al- 
kali doped picene^ Even though the interplay between 
phonon-mediation and electronic correlations is still be- 
ing discussed to better rationalize the superconducting 
transition all across the fullerides family^ the magni- 
tude of the electron-phonon coupling (EPC) in Cqo has 
been the subject of numerous theoretical and experimen- 
tal studies since the early 90s in order to evaluate in par- 
ticular the effective phonon-mcdiated attractive potential 
V ep central to the BCS theory. 

Of particular relevance for electron-doped fullerenes, 
the coupling to the lowest unoccupied molecular or- 
bital (LUMO) was explored extensively on the basis 
of various theoretical approaches^ from earlier combi- 
nations of semi-empirical and density functional the- 
ory (DFT) calculations^— to fully first-principles DFT 
studies*^ - — Concerning the contribution of the H g vi- 
brational modes, values from 38 meV to 68 meV were 
calculated within DFT and (semi) local functionals such 
as the local density approximation (LDA), while the 
A g modes were found to provide a much smaller con- 
tribution, consistently below about 10 meV. These cal- 
culated energies are significantly lower than the avail- 
able experimental values extracted from photoemission 
(PES) experiments on isolated fullerenes with V ep found 
to extend from 96 meV to 147 meV for the H g modes 
contribution; 2 ^— and from 107 meV to 158 meV includ- 
ing both A g and H g contributions! 23 ' 24 

In recent worker— the EPC in 6*60 was revisited us- 
ing DFT and hybrid functionals. An important outcome 



of these studies was a significant increase of V ep with 
increasing percentage of exact exchange within modified 
B3LYP or PBE functionals. This clearly indicates that 
in such systems, not only the electronic excitation en- 
ergies, but also the EPC constants, are very sensitive 
to the choice of the exchange-correlation functional. De- 
spite the overall better agreement with experiments when 
hybrid functionals are used, it is unclear which amount 
of exact exchange should be used for the fullerenes, or 
any finite or extended system in general. Further, the 
evaluation of the coupling constant to individual energy 
levels such as the t\ u state in Cqo relies on the identifica- 
tion of the Kohn-Sham eigenstates and eigenvalues with 
proper electronic quasiparticle states and excitation en- 
ergies. From a pragmatic point of view, the amount of 
exact exchange in e.g. B3LYP is adjusted to reproduce 
ground-state properties of a set of molecular systems^ 
but does not guaranty that Kohn-Sham eigenvalues re- 
produce correctly quasiparticle energies. For example, 
the C 60 Kohn-Sham HOMO-LUMO gap is 2.8 cV within 
the DFT-B3LYP approach^ 2 ! This is better than the 
1.6 eV obtained within DFT-PBE^ but still significantly 
smaller than the 4.9 eV experimental gap in the gas 
phased 

In the present work, we study the electron-phonon cou- 
pling in the Cqo fullerene using the first-principles GW 
approach providing well-defined and accurate quasipar- 
ticle energies within a parameter-free many-body per- 
turbation theory framework. We focus on the threefold 
t\ u lowest unoccupied molecular orbital (LUMO) which 
forms the conducting states in electron-doped fullerides 
and thus determines the superconducting properties. We 
find that the electron-phonon potential V ep is increased 
by as much as 48% as compared to DFT-LDA calcula- 
tions, bridging the gap with experimental data. In par- 
ticular, the contribution from the H g modes is now found 
to be within 4% of the two most recent experimental esti- 
mates. The present results may invite to reconsider pre- 
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vious DFT calculations of the electron-phonon coupling 
constants involved e.g. in the study of the superconduc- 
tivity in molecular or extended systems. 

II. METHODOLOGY 

In the GW quasiparticle formalism,— ~— for which 
decades of expertise exist in the case of bulk systems , 34 ' 35 
the exchange-correlation potential is described by a non- 
local energy-dependent self-energy T,(r,r'\E) which can 
be expressed as follows: 

£ GM/ (r,r'|£) = — [ (L)G(r,r'\E + u))W(T,r'\u)) 
2tv J 

G(r,r'\w) = Y,M*Wn{Y')l(u-e n ±i5) 

n 

W(r,r'\uj) = Jdr" V c (r, r'^e" 1 ^", r'|w) 

where G is the time-ordered Green's function^ and W 
the dynamically screened Coulomb potential built from 
the bare Coulomb potential V c and the non-local inverse 
dielectric matrix e _1 at finite frequency. For extended 
solids, the "starting" (e„,0„) eigenstates used to build 
G and the dielectric response are traditionally obtained 
from a ground-state DFT calculation with (scmi)local 
functionals such as LDA or PBE. 

In the case of isolated molecules, the GW approach 
was recently thoroughly validated on a large set of small 
molecules^ and larger organic ones such as fullerenes, 
porphyrins^ 8 - or DNA/RNA nuclcobases.— An excellent 
agreement with experiment for the ionization energies 
and electronic affinities were obtained through a simple 
self-consistency on the eigenvalues with DFT-LDA eigen- 
states used as the starting point £2t£2. In this latter ap- 
proach, labeled GW in what follows, the HOMO-LUMO 
gap of gas phase Cqq was calculated to be 4.91 eV;2& in 
much better agreement with experiment than the DFT- 
B3LYP Kohn-Sham value. 

Our calculations are based on a recently develope d 38 i 41 
gaussian-basis implementation of the GW formalism 
(the Fiesta code) with explicit treatment of dynam- 
ical correlations through contour deformation tech- 
niques. We start from DFT-LDA eigenstates calcu- 
lated with the Siesta package^ 2 - and a double-zeta 
plus polarization (DZP) basis^ 3 - for the description of 
the valence orbitals combined with standard norm- 
conserving pseudopotentialsi^ 4 - As shown below, the re- 
sulting electron-phonon coupling potentials are very sim- 
ilar to that obtained with all-electron calculations ; 16 ' 19 ' 20 
at least at the DFT level for which several studies 
are available. While GW calculations exploiting DFT 
eigenstates generated with pseudopotentials represent 
the most common approach 3 ^— , a specific aspect of the 
present gaussian-basis implementation is that the auxil- 
iary basis described here below has been optimize d 38 ' 39 ' 45 



to project onto the products of occupied/unoccupied 
pseudized eigenstates. 

The needed two-point operators such as the dynami- 
cal and non-local free-electron susceptibility xo( r , r '|w), 
the screened Coulomb potential W(r, r'|u;) and the 
self-energy operator T, GW (r, r'|o;), are expressed on 
an auxiliary even-tcmperered gaussian basis consist- 
ing of 4 gaussians per each (s,p,c/)-channel with lo- 
calization decay constant (0.2,0.5,1.25,3.2) a.u. Such 
a basis was thoroughly tested in the GW study of 
a retinal chromophore^ of fullerenes, porphyrins or 
phtalocyanincs 3 - and DNA/RNA nuclcobascsi 39 ' 46 

For numerical accuracy when calculating the corre- 
lation contribution to the self-energy, we first evaluate 
T,f w (E) on a coarse energy grid to get a first estimate 
of the quasiparticle energy, and then recalculate T l f w (E) 
on a fine energy grid around this energy to refine our cal- 
culated correlation contribution. We verify as well that 
performing the imaginary-axis integration needed in the 
contour deformation technique (see Ref.HH) with 12 gaus- 
sian points yields results within 0.1 meV as compared to 
a calculation using 20 gaussian points. 

Following the results of Ref. we use the relaxed 
structure and phonon eigenmodes generated within the 
DFT-B3LYP approach and a 6-311G(d) basis4i This ap- 
proach was shown to yield the best eigenfrequencies as 
compared to Raman experiments! 48 ' 49 The EPC matrix 
elements are evaluated using a direct frozen-phonon tech- 
nique. Namely, we deform the molecule along the e v vi- 
brational eigenvectors with typical amplitudes of 0.05 A 
and compute the slope (e„ • V)e, of the variation with 
respect to the deformation amplitude of the DFT Kohn- 
Sham eigenvalues, and further of the GW quasiparticle 
energies. We verify that we remain in the linear regime 
as confirmed by the small value of the error on the re- 
gression coefficient within the fitting procedure. Group 
theory analysis shows that the (t\ u ® t\ u ) direct product 
projects only onto the non-degenerate A g modes and the 
five- fold H g vibrational modes, significantly reducing the 
number of matrix elements to be calculated. It remains 
that ten modes can contribute to the coupling, so that a 
very large number of GW calculations are needed. 

To conclude this methodology section, we note that 
in the present approach based on frozen-phonon tech- 
niques where atoms are explicitly displaced, the calcu- 
lated electron-phonon coupling potentials may be sub- 
ject to errors related to the use of localized basis ( "Pulay 
errors"). This issue was previously explored at the DFT 
level by comparing localizcd-basis and plancwave calcula- 
tions showing small differences (see Supplementary Ma- 
terials, Ref.l20l). In the present case of GW calculations, 
we verify here below that increasing both the size of the 
DFT and auxiliary basis, and taking more diffuse aux- 
iliary orbitals, hardly changes the value of the coupling 
constants, suggesting again small errors related to the 
use of finite atom-centered bases. 
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III. RESULTS 

Our results are gathered in Table I where we provide 
the EPC potential contribution V£ p for each of the ten 
relevant modes, including their degeneracy, namely: 

v e P = _9^_ V l<^l(g.-V)^ C ' F |0j>| 2 

where (e„ ■ \7)V SCF is the normalized variation of the 
self-consistent potential under distortion of the molecule 
along the vibrational mode with index (v), degeneracy 
g v and frequency u„. The indices run over the t\ u 
manifold with <7ti«=3 degeneracy. The above formula is 
the molecular limit^^ of the central definition used in 
ah initio studies of phonon-mediated superconductivity 
in extended solids. 

In the present frozen-phonon approach.— ~— the ex- 
plicit deformation of the molecule diagonalizcs the eigen- 
states with respect to the perturbation, leaving only 
the intraband transitions which, thanks to Hellmann- 
Feynman theorem, can be calculated through the vari- 
ation of the corresponding energy levels, namely: 

3 3 

J2 I < &\(ev ■ V)V SCF \^ > | 2 = J2 lte< ■ V)^| 2 

i— 1 

with derivatives calculated through finite differences. 
This connects EPC matrix elements and the variation 
of the electronic energy levels with respect to vibra- 
tional displacements. This approach is similar to that 
of Refs. Il6l - t20l but we use both the GW quasiparticle 
energies and the DFT Kohn-Sham eigenvalues, allow- 
ing direct comparison. As an internal accuracy test, fol- 
lowing early group symmetry analysis^ the trace of an 
H g perturbation is zero in the (t\ u ) subspace, namely: 
Si=i(eV • V)ffj = 0, a condition which is well verified 
within our DFT and GW calculations. 

Our LDA data yield a total 73.4 meV coupling, in 
good agreement with the 75.8 meV all-electron PBE 
value of Ref.[2(| Comparing to other similar calculations, 
namely extracting the coupling coefficient from the evo- 
lution of the DFT-LDA Kohn-Sham eigenvalues under 
molecular distortion, our 65 meV value for the H g modes 
contribution is also very close to the 68 meV value by 
Antropov and coworkers^ within a full-potential frame- 
work, and the 67 meV obtained with an all-electron gaus- 
sian basis*^ Consistently with early Raman analysis^ all 
studies agree on the predominance of the two high energy 
H g (8) and H g (7) tangential modes, but contributions at 
lower energy such as from the H g (2) and H g (3) radial 
modes are found to be important as well. 

The central result of the present study is the dramatic 
48% increase of the total coupling potential within the 
GW approach as compared to LDA calculations. V ep is 



indeed found to increase from 73.4 meV (LDA) to 108.6 
meV (GW). For the H g modes, the calculated GW value 
of 92.8 meV agrees well with the two most recent 96.2 
meV and 96.5 meV independent experimental estimates 
of Ref. [HI (Table V) and Ref. H, respectively. Further, 
the total GW coupling of 108.6 meV is in close agreement 
with the latest 106.7 meV experimental valuei 21 ' 23 i 50 The 
present results clearly question the accuracy of the EPC 
calculated within DFT and (semi)local functionals. In 
view of the remarkable agreement with experiment ob- 
tained with the present parameter-free GW formalism, 
one can hope that this approach will improve our under- 
standing of phonon-mediated processes in general. 



IV. DISCUSSION 

We can now comment on the recent studies performed 
with hybrids functionals. Since both the experimental 
and GW total coupling potential fall within the rather 
large [93,111] meV energy range obtained by changing 
the amount of exact exchange from 20% to 30% in the 
hybrid DFT approaches ^r— one could certainly build a 
functional yielding excellent agreement with experiment 
for this specific Cgo system. However, clearly, the amount 
of needed exact exchange may vary from one system to 
another (see below). It is further instructing to compare 
mode by mode the various approaches. Considering e.g. 
the A g (2) and H g (8) modes showing large coupling, it ap- 
pears that the largest amount of exact exchange tested 
so far (30%) is not enough to reach the GW results. In 
contrast, the GW values for the H g (n = 2, 3, 5, 6) modes 
are well within the hybrid functionals range. This sug- 
gests that even for a given single molecule, it seems dif- 
ficult to optimize the amount of exact exchange so as to 
reproduce the GW results mode by mode. This observa- 
tion leads to emphasizing the importance of the non-local 
and dynamical correlation part of the GW self-energy. In 
our approach where only the energy levels are updated, 
but not the wavefunctions, the differences between DFT- 
LDA and GW results can only stem from the replace- 
ment of the exchange-correlation functional by the GW 
self-energy. 

An interesting observation is that a non-self-consistent 
Go Wo calculation starting from LDA eigenstates (see col- 
umn 6 of Table I) leads to a coupling constant which is 
still significantly larger than the DFT-LDA value, but 
smaller than the GW one, and very similar to that of 
the hybrid B3LYP functional. As emphasized in recent 
worker— in the case of molecular systems, the signifi- 
cantly too small starting LDA gap leads to a large over- 
screening in the standard Go Wo (LDA) approach. In the 
present Ceo case, the Go Wo (LDA) gap is found to be 
4.4 eV, much better than the 1.6 eV LDA value, but 
still smaller than the 4.9 eV experimental and GW val- 
ues. Qualitatively, this overscreening certainly softens 
the variations of the ionic and electronic potential seen 
by the electrons upon lattice distortion. 
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Theory 






Experiments 




uj (cm i ) 


LDA 


B3LYP 


Hybrids 


Go Wo (LDA) 


GW 


Iwahara a 


Hands" 


Gunnarsson c 


A g (l) 

A a m 

H g (l) 
H g (2) 


496 
1492 
265 
435 


0.5 
7.7 
5.1 
9.9 


1.2 
10.9 
5.8 
10.8 


1.2 - 1.7 
10.5 - 12.9 

5.3 - 6.0 
10.8 - 13.8 


1.2 
13.6 
4.4 
15.3 


I. (107%) 
15.0 (93%) 
6.4 (27%) 

II. 2 (14%) 








H g (3) 
H g (4) 


721 

785 


9.1 
4.2 


11.9 
5.2 


11.0 - 16.7 
4.2 - 5.3 


12.3 
4.7 


13.9 (53%) 
5.6 (36%) 








H g (5) 
H g (6) 


1123 
1265 

1/1/10 


4.2 
2.1 

1 « Q 


5.0 
2.1 

o°. n 


5.0 - 6.7 
2.0 - 4.2 

no n 07 7 


4.2 

2.3 
on n 


5.2 (23%) 

2.3 (9%) 

07 « (f<^7^\ 








H g (8) 


1608 


13.7 


17.7 


17.0 - 19.3 


15.6 


20 .4 (49%) 








Total A g 




8.2 


12.2 


12.1 - 14.3 


14.8 


16.0 (95%) 


10.6 




11.3 


Total H g 




65.2 


81.5 


80.0 - 96.3 


78.8 


92.6 (42%) 


96.2 


96.5 


147.0 


Total 




73.4 


93.7 


93.7 - 110.7 


93.6 


108.6 (48%) 


106.8 




158.3 



TABLE I: Calculated electron-phonon coupling contributions to V ep for the A g and H g modes (meV) calculated within 
LDA, B3LYP (Ref. l2(ih . DFT with various hybrid functionals at the same 6-311G(d) level (column Hybrids with data from 
Refs. Il9l - [2ll ). non-self-consistent GoWo(LDA) and GW. The percentage of increase as compared to LDA is indicated in 
parenthesis. The experimental data are compiled in the three last columns. 



a Ref. 
h Ref. 
c Ref. 



Table V. 



Recently^ the EPC matrix elements in graphene for 
the T-E2 g and K-A[ phonon modes were studied within a 
non-self-consistcnt GoWo(LDA) approach. As compared 
to DFT-LDA calculations, the square of the deforma- 
tion potentials, labeled (-Dp) and (D^), were shown to 
increase by 41% and 114% respectively^ This is consis- 
tent with our own results, suggesting that EPC matrix 
elements are significantly affected by the GW correction 
in both finite and extended systems. A further impor- 
tant outcome of this study was that, in graphene, the 
DFT-B3LYP approach yields significantly too large cou- 
pling constants as compared to experiment, in contrast 
with the present case of fullerenes where the DFT-B3LYP 
calculations underestimate the coupling. This certainly 
points out to the difficulties in obtaining hybrid func- 
tionals which are accurate for both finite and extended 
systems. 



tial V ep increases by 48% as compared to the value cal- 
culated within DFT and (scmi)local functionals such as 
LDA or PBE. The calculated 93 meV GW value for the 
H g modes contribution comes within 4% of the two most 
recent experimental estimates. This demonstrates that 
the present parameter-free approach allows a precise de- 
termination of the electron-phonon coupling potential in 
one of the most studied molecular system. Beyond the 
important case of the fullerenes, the present results call 
for a reinvestigation of previous DFT-based calculations 
of the electron-phonon coupling in organic systems, and 
possibly as well in "covalent superconducting systems" 
such as in particular MgB2 or doped diamond. Simi- 
larly, the important phonon-induced renormalization of 
the electron and hole band width in organic semicon- 
ductors may deserve further inspection beyond previous 
DFT calculations. 



V. CONCLUSION 

In conclusion, we have studied using a first-principles 
GW approach the electron-phonon coupling strength in 
the C*6o fullcrcnc, focusing on the t\ u LUMO state of in- 
terest to the superconducting transition in the fullerides. 
It is found that within GW, the electron-phonon poten- 
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